Disorder-induced Majorana metal in interacting non-Abelian anyon systems 
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We demonstrate that a thermal metal of Majorana fermions forms in a two-dimensional system of interacting 
non-Abelian (Ising) anyons in the presence of moderate disorder. This bulk metallic phase arises in the p = 
quantum Hall state when disorder pins the anyonic quasiparticles. More generally, it naturally occurs for various 
proposed systems supporting Majorana fermion zero modes when disorder induces the random pinning of a 
finite density of vortices. This includes all two-dimensional topological superconductors in so-called symmetry 
class D. A distinct experimental signature of the thermal metal phase is the presence of bulk heat transport down 
to zero temperature. 
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Quantum matter with exotic non-Abelian quasiparticle 
statistics has become a highly sought-after state of matter 
IT]. Theoretical proposals suggest the existence of the most 
elementary incarnation of such a state - arising from zero 
modes of Majorana fermions - in diverse two-dimensional 
(2D) systems such as the fractional quantum Hall liquid at 
V = 5/2 filling fT], ^x+^Vy superconductors |3|, heterostruc- 
tures of topological insulators and superconductors |4|, and 
possibly certain Iridates which effectively realize the Kitaev 
honeycomb model |5|. However, despite intense efforts, ex- 
perimental observation of these states has proven elusive to 
date. In this manuscript we point out a distinct collective 
phenomenon of these states in the presence of disorder - the 
formation of a thermal metal manifesting itself in bulk heat 
transport which may be amenable to a more direct observa- 
tion than previously considered experimental signatures such 
as sophisticated interferometric setups 16 |. As we show, this 
thermal metal phase aiises when a finite density of quasipar- 
ticles is spatially pinned by disorder, and the resulting macro- 
scopic degeneracy - characteristic of non-Abelian statistics - 
hybridizes in the presence of interactions between the quasi- 
particles. For concreteness, we will discuss our results in the 
context of the v = b/2 quantum Hall state, but point out that 
our results hold for other 2D systems with randomly pinned 
zero modes of Majorana fermions. 

The disordered v — b/2 quantum Hall state.- Since the 
seminal work by Moore and Read \1\ the quantum Hall state 
on the = 5/2 plateau has been a prominent candidate state 
supporting quasiparticles which obey non-Abelian statistics. 
Detuning the magnetic field from the middle of this plateau 
will introduce a finite density of charged quasiparticles (or 
quasiholes) which for a clean system will naturally form a tri- 
angular Wigner crystal, as recently observed for the = 1/3 
quantum Hall state |8| . Even though pinned, these quasipar- 
ticles still carry a fluctuating quantum mechanical degree of 
freedom, similar to the spin degree of freedom of an electron 
in a Mott insulator. It is this degree of freedom, which encodes 
the non-Abelianness of the quasiparticle, that is referred to as 



an Ising anyon a. Similar to a pair of ordinary spin- 1/2 mo- 
ments, which can combine as 1/2x1/2 = + 1 into a singlet 
or triplet state, two Ising anyons can form two distinct collec- 
tive states 

cr X a- = 1 + , (1) 

where the state 1 refers to the vacuum (or trivial) state and ijj is 
another quasiparticle excitation. An immediate consequence 
of this so-called fusion rule ([T]) is the formation of a macro- 
scopic degeneracy of states, which in the case of N ordinary 
spin-1/2 moments grows like 2^ and for the case of Ising 

/-N 

anyons grows asymptotically as v 2 for a system of N Ising 
anyons. This macroscopic degeneracy is lifted by interactions 
between the Ising anyons. The most elementary example of 
such a splitting occurs for a pair of Ising anyons brought into 
spatial proximity which will naturally split the energies of the 
two states on the right-hand side (rhs) in Eq. ([T]i, in analogy to 
the splitting of singlet and triplet states for a pair of ordinary 
spin-1/2 moments. For more than two anyons the many-body 
problem aiising from these interactions can be formulated |9| 
in terms of pairwise projectors lijk 1 10| energetically favoring 
one of the two states on the rhs of Eq. ([TJ 

H = -2Y,J]k^jk, (2) 

where we choose a convention that lij^ projects onto the triv- 
ial state 1 of Eq. ([T]i. In analogy with ordinary spin systems 
we refer to ,Jjk > as antiferromagnetic (AFM) interactions 
and Jjk < as ferromagnetic (FM) interactions. 

The microscopies governing the interaction strength Jjh are 
not universal and in particular depend on the underlying phys- 
ical system. For the v = 5/2 quantum Hall state it has been 
shown 1121 that the interactions between quasiholes (qh) ex- 
ponentially decay with the separation between the qhs (on the 
scale of the magnetic length) and oscillate in sign (analogous 
to RKKY interactions between magnetic moments). Simi- 
lar behavior, though quantitatively different in the sign and 
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FIG. 1. Going from Ising anyons (left) to Majoranas (right). 

strength of the oscillations, was found for vortex interactions 
in Paj+ipj, superconductors 1 13| and in the Kitaev honeycomb 
model [14|. 

In the presence of spatial disorder, the Wigner crystal does 
not form an ideal triangular lattice but rather exhibits a ran- 
dom distribution of bond lengths. If the variation of these 
bond lengths exceeds the characteristic length of the sign os- 
cillations (the magnetic length of the quantum Hall state), 
the exchange coupling constants Jj^ will exhibit significant 
sign disorder. Due to the exponential decay of the coupling 
strength with distance, the essential physics of the anyon- 
model with random interactions is expected to be captured by 
a Hamiltonian of the form (|2| with independently distributed 
nearest-neighbor exchange interactions Jjk of random sign. 
Recalling the analogy of the anyonic degrees of freedom to 
ordinary magnetic moments, the random anyon problem thus 
bears a striking similarity with the Edwards-Anderson quan- 
tum spin glass [15|. 

One natural angle to study the random anyon model is to 
pursue a strong-randomness renormalization group (RG) ap- 
proach |16|. However, as we will report elsewhere |17|, the 
system flows towards weaker disorder under this RG and thus 
is not governed by an infinite randomness fixed point. This 
renders the strong-randomness RG approach incapable of cap- 
turing the essential physics of the problem. Below, we there- 
fore follow an alternative path exploiting the equivalence of 
our random anyon model with the problem of Anderson lo- 
calization of disordered Majorana fermions. This problem is 
accessible to exact diagonalization studies of large samples 
which allow us to identify the analytic low energy theory for 
the underlying phase. 

Majorana fermions.- For the v = 5/2 quantum Hall 
state, quasihole (or quasiparticle) excitations are Ising anyons, 
which can also be thought of as vortex-like excitations carry- 
ing edge degrees of freedom as illustrated in Fig.[T| The non- 
Abelian aspects of the Ising anyon <t can be described in terms 
of the Majorana fennion field living on this edge. The Majo- 
rana fermion on an isolated vortex has an exact eigenstate at 
zero energy. We denote the creation operator for such a Ma- 
jorana fermion zero-mode by jj. This represents the Ising 
anyon at position j. In terms of these Majorana operators the 
original anyon Hamiltonian (|2]) can be recast into the form 

= ^ ij]kl]lk (3) 

{jk) 

up to an additive constant. Because Majorana fermions are 
their own Hermitian conjugates, i.e. 7j = we can in- 
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FIG. 2. Phase diagram of sign-disordered triangular lattice hopping 
problem in the (p, ^) plane. 

terpret Hamiltonian (|3]l as a hopping problem where the ma- 
trix Jjk must be real and antisymmetric in order to ensure 
Hermiticity 1 18 1. Thus the interacting anyon problem reduces 
to a hopping problem of non-interacting (Majorana) fermions 
lfT9]| . This mapping was used to solve the uniform interact- 
ing anyon problem on a triangular lattice |20|, which maps 
to a hopping problem with uniform 7r/2 flux around each tri- 
angular plaquette, producing a gapped topological state with 
non-vanishing Chern number. This was later understood to be 
a special case of topological liquid nucleation \2li22\ in inter- 
acting non-Abelian systems. In general, two distinct gapped 
topological phases emerge depending on the sign of the (uni- 
form) coupling J in any given non-Abelian theory 11211 l22l . 
For the case of Majorana fermion zero-modes, the nucleated 
topological liquids are both Abelian and if brought into spa- 
tial proximity would be separated by an edge state of confor- 
mal central charge c — 1, which corresponds to a difference 
Ai/ = 2 in Chern numbers. 

In the presence of disorder we are interested in the physics 
of anyons with interactions of random signs. We introduce 
sign disorder |23| by randomly flipping the sign of the cou- 
plings Jjk = ±1 on the lattice independently with probability 
p G [0, 1]. The limits p ~ and p = 1 then correspond to 
the uniform problem with J = +1 and J = — 1, for which 
the ground states are gapped topological states, which can be 
characterized by Chern numbers ly = —1 and = +1, re- 
spectively. Around these limits these topologically-insulating 
ground states appear to be stable to a small amount of sign 
disorder, as indicated in the phase diagram of Fig. [2] They 
are gapped only in the limit of strictly zero disorder, while 
for weak disorder there should be a low density of Anderson- 
localized states for energies \E\ less than a critical energy. We 
note that the phase diagram is symmetric under p 1 ~ p 
followed by time reversal of the Majorana fermions, which is 
a (statistical) symmetry of the ensemble of Hamiltonians ([3]l. 

In order to determine what happens at intermediate disor- 
der, we numerically investigate the fermion hopping problem 
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FIG. 3. Finite-size scaling collapse of the disorder averaged Chem 
number at ^ = —0.1 with pc ~ 0.225 and « 5 ± 1. 
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FIG. 4. Crossover of the disorder averaged Chem number i/ at /j, = 
with no visible finite-size scaling. 



The thermal metal- In the language of Anderson local- 
ization, the random Majorana hopping problem ([3]) falls into 
what is known as symmetry class D 1251 l26l . It has been 
shown llT9l l27l ^8 1 that problems in this symmetry class can 
exhibit three distinct phases: a topological insulator (or ther- 
mal quantum Hall insulator), a (topologically) trivial insula- 
tor and a thermal metal phase. The most remarkable of these 
three phases, the thermal metal phase, has previously been 
found only within the framework of Chalker-Coddington net- 
work models ||281|29I , which are not directly related to a mi- 
croscopic situation. We will show that the v = Q state oc- 
curring at vanishing /i = in our random Majorana hopping 
problem around p = 1/2 is precisely this thermal metal phase 
1301, thereby also providing the first concrete microscopic sit- 
uation for the appearance of this metallic state. 

The most direct evidence for this metallic state comes from 
the density of eigenstates of Hamiltonian ([3]) near zero en- 
ergy E — Q obtained in the numerical solution. We have 
calculated disorder averages of the density of states p{E) 
using standard sparse diagonalization techniques for system 
sizes up to 256 x 256, again with ^ 2 x lO"* disorder real- 
izations for a given value of disorder strength p. These are 
plotted on a semilog scale in Fig. |5] for p = 0.5, clearly 
revealing a logarithmic divergence characteristic of the pre- 
dicted behavior for the thermal metal phase in class D f32l. 
Further, the bump and oscillations visible at the lowest en- 
ergies (of the order of the inverse mean level spacing) di- 
rectly follow the form predicted from random matrix theory 
125 J p{E) — a + sin{2'!TaEL'^)/{2TrEL'^) , where a is anon- 



(|3]l. To this end, we diagonalize the matrix ij^jk from Eq. (|3]l 
by exact diagonalization techniques for lattices ranging from 
L X L = 8 X 8 to 64 X 64 and considering 10'^ — 10** disorder 
realizations for various disorder strengths p. For each disorder 
realization we calculate |24| the Chern invariant h' for eigen- 
states filled up to a chemical potential /i. While it is benefi- 
cial to consider the phase diagram in this enlarged parameter 
space {p, /i), we note that the physically relevant ground state 
of Hamiltonian (|3]l corresponds to the line of precisely /i = 0. 
Going to non-zero /i first, we find a clear transition from a 
topological insulator with ly = —1 to a topologically trivial 
phase with = as p is increased from to 1/2. We can esti- 
mate the critical coupling Pc(m) for this transition via a finite- 
size scaling collapse of the disorder averaged Chern number 
P, as shown in Fig. [3] These estimates provide the indicated 
phase boundary of the topological insulator in the phase di- 
agram of Fig. |2] At /i = 0, the chemical potential relevant 
to the original Majorana problem ([3]), the disorder averaged 
Chern number P shows a slow crossover from = — 1 in the 
topological insulator to a state with v = Q around p = 1/2 
with no visible finite-size scaling, as shown in Fig. |4] This 
state with j/ = at vanishing chemical potential /i = is not 
the trivial insulator seen at finite chemical potential, but rather 
a metallic state as we will demonstrate in the following. 
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FIG. 5. Disorder averaged density of states p{E) near iJ = at 
p — 1/2 collected from = 2 — 3 x 10^ samples per size L. 
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FIG. 6. Left: The inverse participation ratios of the wavefunction 
for disorder strength p = 0.5 averaged over N — 2 x 10'^ samples 
per size L in an energy window E < 10^'^. Right: The multifractal 
exponents characterizing the decay of the inverse participation ratios. 

universal parameter. 

Further evidence of the metallic state reveals itself upon 
close inspection of the wavefunctions near E = 0. A metal- 
lic state in symmetry class D has a weak antilocalization cor- 
rection corresponding to a marginally irrelevant perturbation 
|32|, in the RG sense. This implies that the moments of the 
wavefunctions, or the so-called inverse participation ratios, 
exhibit 'weak multifractal' behavior (for the q-th moment) 

with multifractal exponent Tq = 2{q — 1) — ''yq{q — 1), where 
the non-universal coefficient 7 depends only logarithmically 
on the linear system size L. The moments of our numerically 
obtained wavefunctions precisely follow this predicted decay 
with system size as shown in the left panel of Fig. |6] The 
extracted multifractal exponents for different moments q are 
fitted by = -O.lSlg^ + 2.l35q - 1.979, which matches 
the predicted one-parameter form with 7 = 0.151. Finally, 
we note that the thermal metal is robust, in particular with 
respect to potential higher-order interactions between the Ma- 
jorana zero modes ji in Eq. (|3]l. These have been found to be 
irrelevant (in the RG sense) in the metallic phase [Mi- 
Conclusions.- To summarize our results, we have demon- 
strated that moderate disorder induces a thermal metal phase 
in systems which harbor a finite density of interacting Majo- 
rana fermion zero modes. A prominent candidate in which to 
look for this thermal metal phase may be the i/ = 5/2 quan- 
tum Hall plateau, which has long been suggested to exhibit 
non-Abelian anyon excitations. In fact, the experimental ob- 
servation of thermal heat transport in the bulk may help to re- 
veal the existence of this non-Abelian state. In abelian states, 
the thermal conductivity (divided by temperature) would van- 
ish when approaching zero temperature, in contrast to the log- 
arithmic divergence exhibited by the disordered thermal metal 
II32I . The occurrence of this thermal metal state is not re- 
stricted to the v ^ 5/2 quantum Hall state, but is expected 
to be present in any two-dimensional time reversal breaking 



topological phase supporting pinned Majorana zero modes 
in the presence of moderate disorder - this includes, e.g., 
Px+iPy superconductors and microscopic realizations of Ki- 
taev's honeycomb model. 
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